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Wave Functions for Open Quantum Systems 
and Stochastic Schrodinger Equations 
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It is shown that evolution of an open quantum system can be exactly described in terms of wave 
function which obeys Schrodinger equation with randomly varying parameters whose statistics is 
universally determined by separate dynamics of the system's environment. Corresponding stochastic 
evolution of the wave function is unitary on average, and this property implies optical theorem for 
inelastic scattering as demonstrated by the example of one-dimensional conducting channel with 
thermally fluctuating potential perturbation. 
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I. Generally, even simple quantum systems have no 
definite wave function and require the density matrix 
language instead 0. But if the density matrix, R{t) , 
evolves under von Neumann equation, still one can write 



(1) 



with a{t) obeying the Schrodinger equation, and say 
that the system has random wave function whose ran- 
domness is determined by initial conditions only. This is 
not the case as soon as the system becomes open, i.e. con- 
stantly interacts with the rest of the nature, and there- 
fore is governed by some kinetic equation. However, in 
the framework of the "stochastic representation" of quan- 
tum interactions or deterministic interactions at 
all ^ 1^ , we can again rehabilitate wave function if say 
that it is stochastic one being governed by a well-defined 
stochastic Schrodinger equation. Let us consider this po- 
tential convenience more attentively. 

II. Suppose that (i) a composite quantum system 
"D-l-B" , with "D" being a dynamic subsystem under di- 
rect interest and "B" its environment ( "thermal bath" or 
'all the Universe" or some other), has the bilinear Hamil- 
tonian 



H{t) - Hd{t)+Hb 



(2) 



(repeated indices imply summation over them), where 
operators Ha , Dj relate to "D" while Hb , Bj to "B" ; 

(ii) common density matrix of "D-l-B" , p{t) , obeys the 
von Neumann equation, dp(t)/dt = i[p(t) , H (t)] / h , and 

(iii) in the old days "D" was statistically independent 
of "B", p{to) = PD{to) X psito) (where e.g. ^ -oo ). 

Then marginal density matrix of the subsystem "D" , 
PD{t) = Ttb p{t) , can be represented 0,01 as statistical 
average, 



pnit) - {R{t)) , 



(3) 



of a random density matrix, R(t) , which satisfies the 
stochastic von Neumann equation 



dR 
'dt 



= -{[R,Ha{t)]+^*{i)RDj-mDjR] , (4) 



where S,j{t) are complex random processes which com- 
mute one with another and with any other objects^ the 
superscript * means complex conjugation and (...) sta- 
tistical average with respect to S,j{t). Sometimes it is 
convenient to write 

e(t) = x{t) + iny{t)/2 , r(<) = x{t) - ihy{t)/2 , 

with Xj (t) and yj (t) treated as real- valued random pro- 
cesses. Statistics of ^{t)^s is completely determined by 
separate dynamics of the subsystem "B" : 



( exp{ / [ a* {t)C, it) + a, m* (t) ]dt}) = 



(5) 



Tr^ ex]3 



ajit)B,{t)dt 



5xp 



a*At)Bj{t)dt 



PbQ 



with Bj{t) = exp[iHi,t/h] Bj exp[—iHi,t/h] , left ar- 
row and right arrow denoting chronological and anti- 
chronological time ordering, respectively, aj{t) and 
a*{t) being arbitrary test functions, or probe functions 
(which are not necessarily mutually conjugate), and 
Pbo = exp[iHbtQ/h] psito) exp[-iHbto/h] . 

Clearly, any solution to Eq^can be represented in the 
form (P), where wave functions (quantum states) ^'^(i) 
satisfy the stochastic Schrodinger equation 



--{Ha{t)+Ut)D,+<t)}-^{t) (6) 



with initial condition ^'0,(^0) ^'q ^l(io) = Poito) , and 
e{t) is a suitable (real-valued) gauge function. At that, 
formula (jsj supplies exhaustive statistical information 
about the noises ^j{t). To avoid a need for e{t) , of 
course, it is reasonable to assume 
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Ti-D Dj=Q , Tib Bj=0, 



(7) 
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thus excluding from ^ trivial contributions without fac- 
tual interaction. 

III. Because of the complexity of £,j{t), the stochastic 
evolution under EqElis not unitary: 



- {^^it)\^pit)) = y,it){^^{t)\D,\^pit)) + (8) 
But this evolution is unitary on average : 



dt 



{^a.{mp{t)) ) = {y,{t){^4t)\D,\^p{t)) ) = , (9) 
( i^o^impit)) ) = (*„(io)|«'/3(io)) « So^p (10) 



The statement follows from Q and 'i>a{to)'s defi- 
nition, while the statement ^ from remarkably singular 
statistical properties of yj{t) 

To see the ?/j(i)'s singularity, it is sufficient to as- 
sign in ISJ a{t) — iu{t)/h and a*{t) ~ —iu{t)/h when 
t > 9, with real test functions u{t) and 6 > to. 

Then ^ [J a* {t)B{t)dt] = U ^ f a*{t)B{t)dt and 



ex]5 [/ a{t)B{t)dt] 



e5q5 



f a{t)B{t)dt W , where 



U ~ fexp [~i fgu{t)B{t)dt/h\ . Since U is unitary, 
U'^U = 1 , in fact it disappears from r.h.s. of (O: 

(exp{ /, u{t) y{t) dt} exp{ f[a*mit) + a(t)r(t) ] dt}) 
= TrBex^ fa{t)B{t)dt fecp f a*{t)B{t)dt pw 

(11) 

Here from it follows that any statistical correlation be- 
tween y{t) and any earlier variables is zero. In partic- 
ular, the equality lO is valid because 5'(t)'s undergo 
the causality principle and depend on earlier ^(< t)'s 
only. The EqEl implies also {y{ti)...y{tn)) = , 
{y{t)x{t' < t)) = 0, etc. 

At the same time, generally {x{t)y{t' < t)) ^ and 
other y(t)'s correlations with later variables also are 
nonzero. 

Let the noise is stationary, that is [Hi,, pbo] = , and 
besides has zero mean value, {^j{t)) = 0. The latter is 
ensured by the natural condition Tr^ Bj pi,Q = . Then 
0, according to EqI3 



(e*(r)6n(0)) =i^,™(T) , 
K,^{t) = Ti-B B,{T)B,n{0) Pbo = K*^j{-t) 
(e,(r)^™(0)) = (f*(r)C(0))* = 



(12) 



where 9{t) is Heaviside step function. Consequently, 

^ {x,{T)x,n{0)) = ReK,„,{T) , 
K^y^ir) ^ {x,{T)y^m = i2/h)e{T)lmK,„,{T) ^^^^ 

If pbQ is canonical distribution, pi,Q oc e:sq){—Hi,/T) , 
then the identity pi,QB{t)p^^ = B{t -f ih/T) takes place 
and produces additional symmetry relation 



The latter if combined with H12() results in the 
"fluctuation-dissipation relations" between spectral com- 
ponents of correlators if^^ (r) and Kj^ (r) as expressed 



by their Fourier expansions: 



KjZi'r) = /o°° cos(wt) aj^{uj)duj , 



K^y (r) = - 



sin(a;T) tanh 



huj 



2T 



(15) 



{uj)duj , 



(16) 

with ajm{^) being some non-negatively defined spec- 
trum matrix. Unfortunately, in |2| formula (|16l) was 
written with wrong plus sign on its r.h.s., and then this 
mistake migrated to |^^|^. But, fortunately, it could 
not influence results of applications of (|16|l since all that 
were even functions of K^]L (r) . The term "fluctuation- 
dissipation" reflects that Kj^ (t) and other cross correla- 
tions between x(t)'s and y{t)'s describe also response of 
"B" to its perturbation by "D" and thus possible energy 
outflow from "D" into "B" . For details and extensions of 
the stochastic representation ©-(Ell see 

IV. For illustration, let "D" be one-dimensional con- 
duction channel with finite energy band, that is formed 
by discrete sites s = ... — 1,0,1,..., and its perturba- 
tion by "B" acts as randomly varying potential localized 
at the only site s = 0. Thus ^j{t)Dj = £,{t)D where 
Dss' = 1 if s = s' = and Dss' — otherwise. Defini- 
tion of the own Hamiltonian of "D" is convenient in the 
momentum space: {Hd{t)}kk' = 5{k — k')E{k) , where 
— TT < k < TT and, for instance, E{k) — AE{1 — cos k)/2 
with AE being channel's bandwidth . 

In this example, it is natural to treat ^'(t) in the co- 
ordinate space, choose e{t) = , and transform Eq^linto 
integral equation: 

•i>{t,s) = ^°{t,s) + ^^{t,s) , 
^^{t,s)^-fJlG{t-t\s)m^{t',0)dt' 



(17) 



where 4'°(t, s) and '^^{t,s) are free and scattered parts 
of the wave function, respectively, and 



G(r, s) = /_^^ exp [iks - iE{k)T/h] dk/2n 



(18) 



is free Green function. 

Consider stochastic scattering of a single particle. Let 
initially, at time to = ^^o < , it had wave number 
fco > and velocity vq = dE(ko)/hdkQ, and its wave 
packet took place —L < s < with L « uoTo >> 2n/ko . 
Then we observe the particle at time t > . At that, we 
are interested in the case when duration of the scattering 
process, ~ L/vq , essentially exceeds correlation time, 
Tc , of the scatterer fluctuations ^{t). This is possible if 
L is sufficiently large while vq sufficiently small, at least 
if L > Lc = y/irAETc/h, (for example, if AE = 1 eV 
and each site occupies 3 • 10~* cm then Lc ^ ^/fc where 
Tc is expressed in seconds and Lc in centimeters). 

Obviously, now the equality H10|) can be written as 



K 



j™(t - th/2T) = K„rj{-T - ih/2T) 



(14) 



(EI*(i,5)P)=EI*(io,s)p = 

= EI*°(io,s)P=EI«'°(^,s)p 



(19) 
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Combining this with the first of EqsElone easy obtains 

2[*"*(*5>+*"(*^^*)] + (^|*^^|2) = , (20) 

where all functions are related to the observation time t . 
This formula is nothing but "f/ie averaged optical theorem 
for inelastic scattering'^ (clearly, fluctuations ^(t) in the 
course of scattering can make it inelastic). 

As it was agreed, at time t > the scattering event 
finishes, therefore, we can express final quantum proba- 
bilities of the particle's reflection, 3? , and transmission, 
T , simply as 



5i-E.<ol^l' =E.<ol*^l' -sEI^^P 



(21) 



We took into account also mirror symmetry of the scat- 
tering: the Green function G{t,s) and thus \E''^(i, s) are 
even functions of s (at least when E(k) is an even func- 
tion). Now, taking average of (|21|) and applying (|19|l and 
(Pn|l . we obtain 

{n)=-Rcj:^^*{t,s){^^{t,s)) = l-{T) (22) 

Hence, the optical theorem helps to find the mean re- 
flection probability ( 3? ) directly from mean value of 
\l/(t, s) , without knowing its second-order statistics. 

Another useful expression for ( 3? ) arises if com- 
bine with ((T7jl and then apply the group property 
^ G{ti,s - s')G{t2,s') = G(ti +T2,s). This yields 

{^)=Rej-Jl ^^*it',0){mm)dt' , (23) 

where tp{t) = ^(^,0) satisfies closed integral equation 

m = ^"{t, 0) - i ji Goit - t')at'm') df (24) 

with Go(r) = G(t,0) . 

Alternatively, at i > we can consider amplitudes of 
outgoing waves, 

A,^Al + A^, , Al^S{k~ko) , 

and the probabilities of corresponding particular variants 
of scattering. 



(Pk 



Pk = \Ak\^ 



\{A,)\' + {A,,Al) 



(25) 



Here and below {a , b) = {a b) — (a) {h) is Malakhov's 
cumulant bracket. If ^'''(t, s) was normalized to unit, 
then the EqsEl^re equivalent to 



exp[iE{k){t' - to)/n] £.{t'yj{t') dt' 



(26) 



h^/Z J to 

In terms of A^ , the optical theorem reads 

Re«)=Re(A^,J = -(5R) , (27) 
where again the symmetry is taken into account. 



Let us divide 3? and T into elastic and inelastic parts 
marked by superscripts "el" and "in" . Obviously, 

^ , (T'^O ^ (PfeJ (28) 

These identities together with Eg 1251 and Eal27lvield 

(29) 



(3?™) = (3?) - (3fi'=0 < (3?) - (3fi)2 
(T™) = (T) - {T"') < (T) - (r)2 



In fact, that are one and the same relation, because 
sjjm _ q-m -j-g ^Yie rnirror symmetry. Here from we 
obtain the restriction of total inelastic scattering: 



(3?*") -I- {T'") < 2 (3?) [1 



< 1/2 (30) 



This is consequence of the "unitarity on average" . 

V. Now consider the stochastic scattering more con- 
cretely, basing on the Eg 1241 and assuming that ^(t) is a 
stationary random process, in particular, its mean value 
is time- independent: {(,{t)) — {x{t)) — x = const (of 
course, x must be real). 

Let us use the star * as symbol of causal convolution: 
{a* b){t) = /(^j/ a{t — t')b{t') dt' . Due to the stationarity, 
when averaging Eg 1241 we can write 

mm) = il s(t - t')m'))dt' ^ (s * m{t) , 

= [l + (*/n)Go* S*]-i*0(i,0) 

(31) 

The kernel ^(r) of the "mass operator" looks as 

E{T)^xS{T)~f-G,iTmr),m) + - , (32) 

where dots substitute higher-order "non-factorable Feyn- 
man diagrams" , and 

G,{T) = [l + {u:/n)Go*]-'Go{T) (33) 

Correspondingly, Eql^Hl together with H31|l - H32|) yields 

ReG/(H-G), (T) =Re 1/(1-HG) 



G = (i/hvo) /q°° ^{t) exp{iEoT/h) dr 



(34) 



We noticed that Go{t) expiiEr/h) dr = 1/V{E) , 
where V{E) = y^E{AE — E)/fi is particle's velocity as 
function of its energy (thus uq — y{Eo) )■ 

Because of the stationarity, both {i'it)) and 
{i{t)i}j{t)) oscillate like the source *°(i,0) in it^ : 
(V-'(i)) oc {£_{t)^/j{t)) cx exp{—iEot/'h) . Hence, any in- 
elastic amplitude of H2()|) is zero on average: 



(Af) (X S{\k\ - Ifcol) 



(35) 



From this point, concentrate on simplest but principal 
situation characterized by three assumptions as follow. 

(i) C(^)'s fluctuations are produced by equilibrium 
thermal bath, therefore, according to l|12l) - (|16|l . 

(e(r),C(0)) ^Ki\r\) , {Cir),m) ^ K{t) , (36) 
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K{t) = 



exjpiTiLj /T)e 



1 + exp(ncj/r) 



■ cf{uj) duj , (37) 



where T is bath temperature and <j{lo) > . The 
spectrum <j{lo) can be naturally introduced through 
correlation of the random potential x{t) — Re^(t): 
{x{t),x{0)) = cos{uJT)a{uj)duj. 

(ii) C(t)'s fluctuations are sufficiently small and short- 
correlated (wide-band), e.g. in the sense of 



max cr(w) ^ hAE 



(38) 



so that the scattering can by analyzed in the "Born ap- 
proximation" . 

(iii) X = , i.e. "there is no potential irregularity and 
no scattering on average" . 

Under these conditions, when averaging and other 
variables one can manage with the correlation function 
K{t) only, as if £,{t) was Gaussian noise. Due to (|38|l . 
formula H34|) strongly simplifies to 



1 



Voh J-Eo/h 



V{Eo + fiu) 



Re C « I V./^n . \ , 

(39) 

where F(X) = + cxp(X)] is the Fermi distribution 
function, and presumably Re C < 1 . The contributions 
to C from cj > and uj < correspond to scatter- 
ing with radiation or absorption of energy by the bath, 
respectively. 

When X ^ , then clastic backward scattering must 
disappear at ct — s- , that is {A^i^^) cx u at small a. 
What is for fiuctuational contribution, {A^ko : ^Ifco) ' 
mean probability of this scattering, it is small to the ex- 
tent of (j/L , as well as {Aj, , A^) at all (see Eal40lbelow'). 
Hence, (3?*^') ^ (3?) , and the quantity (|39l) in fact rep- 
resents the probability of inelastic reflection (3?™) . 

The latter statement can be confirmed by direct calcu- 
lation of {Ak , A'^) using Eg 1261 and the approximation 



which is valid under above assumptions. The result is 



{A, , Al) = <7{\s{k)\/n) F{e{k)/T) , (40) 



where e{k) = E{k) — Eq . With use of the equality (jS^Jl, 
the inelastic refiection probability is 



■fco < 



oiAk,Al) 



(41) 



In fact, of course, for the finite incoming wave packet, this 
sum contains a finite number « L only of distinguishable 
outgoing states. The factor 27r/L in (|40|) is just their 
separation in the fc-space, and because of L ^ 1 the 
sum (HU coincides with r.h.s of that is (3?™) « (3?). 

VI. We illustrated once again that the general 
stochastic representation approach to open systems as 
suggested in '2| and developed in can be useful 

in concrete considerations replacing complicated dynam- 
ical evolution of composite systems by stochastic evolu- 
tion of a simple partial system. In the above example, 
if the energy of incoming particle is not fixed but has 
thermally equilibrium distribution with the same tem- 
perature as the scatterer's temperature, then it is easy to 
verify that the same distribution will be reproduced after 
the stochastic inelastic scattering. Hence, we also obtain 
additional evidences that the stochastic representation 
ensures automatical agreement with thermodynamics. 

To conclude, notice that our approach allows to set 
the problems about transport processes and their noise 
(including fiuctuations in transport rates, e.g. in electri- 
cal conductivity), as well as about dephasing of coherent 
quantum processes, in a new fashion, with attraction of 
reach results on linear dynamic systems with randomly 
varying parameters (in particular, on wave propagation 
in random media). 
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